Nonmonotonic effects of parallel sidewalls on Casimir forces between cylinders 
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We analyze the Casimir force between two parallel infinite metal cylinders, with nearby metal 
plates (sidewalls), using complementary methods for mutual confirmation. The attractive force 
between cylinders is shown to have a nonmonotonic dependence on the separation to the plates. 
This intrinsically multi-body phenomenon, which occurs with either one or two sidewalls (general- 
izing an earlier result for squares between two sidewalls) , does not follow from any simple two-body 
force description. We can, however, explain the nonmonotonicity by considering the screening (en- 
hancement) of the interactions by the fiuctuating charges (currents) on the two cylinders, and their 
images on the nearby plate(s). Furthermore, we show that this effect also implies a nonmonotonic 
dependence of the cylinder-plate force on the cylinder-cylinder separation. 

PACS numbers: 12.20.Ds, 42.50.Ct, 42.50.Lc 
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Casimir forces arise from quantum vacuum fluctua- 
tions, and have been the subject of considerable theo- 
retical and experimental interest [U, 0, S 13] ■ We consider 
here the force between metallic cylinders with one or two 
parallel metal sidewalls (Fig. [1]) using two independent 
exact computational methods, and find an unusual non- 
monotonic dependence of the force on the sidewall sepa- 
ration. These nonmonotonic effects cannot be predicted 
by commonly used two-body Casimir-force estimates, 
such as the proximity-force approximation (PFA) 
that is based on the parallel-plate limit, or by addition 
of Casimir-Polder 'atomic' interactions (CPI) [1, 0, S|- 

In previous work, we demonstrated a similar nonmono- 
tonic force between two metal squares in proximity to 
two parallel metal sidewalls, for either perfect or realistic 
metals This work, with perfect-metal cylinders p^ . 
demonstrates that the effect is not limited to squares (i.e., 
it does not arise from sharp corners or parallel flat sur- 
faces), nor does it require two sidewalls. The nonmono- 
tonicity stems from a competition between forces from 
transverse electric (TE) and transverse magnetic (TM) 
field polarizations: In the latter case, the interaction be- 
tween fluctuating charges on the cylinders is screened by 
opposing image charges, in the former case it is enhanced 
by analogous fluctuating image currents. Furthermore, 
we show that a related nonmonotonic variation arises 
for the force between the cylinders and a sidewall as a 
function of separation between the cylinders, a geometry 
potentially amenable to experiment. 

Casimir forces are not two-body interactions: quantum 
fluctuations in one object induce fluctuations throughout 
the system which in turn act back on the first object. 
However, both the PFA and CPI view Casimir forces 
as the result of attractive two-body ( "pairwise" ) interac- 
tions. They are reasonable approximations only in cer- 
tain limits (e.g., low curvature for PFA), and can fail 



qualitatively as well as quantitatively otherwise. Pair- 
wise estimates fail to account for two important aspects 
of the Casimir forces in the geometry wc consider [11 



First, a monotonic pairwise attractive force clearly can- 
not give rise to the nonmonotonic effect of the sidewalls. 
Second, the application of either method here would in- 
clude two contributions to the force on each cylinder: 
attraction to the other cylinder and attraction to the 
sidewall(s). If the latter contribution is restricted to the 
portion of the sidewall(s) where the other cylinder does 
not interpose ("line of sight" interactions), the cylin- 
der will experience greater attraction to the opposite 
side, thereby reducing the net attractive force between 
the cylinders [llj. In contrast, exact calculations pre- 
dict a nonmonotonic force that is larger in the limit of 
close sidewalls than for no sidewalls. These important 
failures illustrate the need for caution when applying 
uncontrolled approximations to new geometries even at 
the qualitative level. (On the other hand, a ray-optics 
approximation, which incorporates non-additive many- 
body effects, at least qualitatively predicts these features 
for the case of two-square/sidewall [Tlj.) 

Figure [T] is a schematic of the three-dimensional ge- 
ometries that wc consider: two infinite, parallel, perfect- 
metal cylinders of radius R separated by a distance a 
(center-to-center separation 2R + a) and oriented along 
the z axis, with one (top panel) or two (bottom panel) 
infinite metal sidcwall(s) parallel to the cylinders and 
separated from both by equal distance h {h + R to the 
cylinder axes). For perfectly conducting objects with z- 
translational symmetry, the electric (E) and magnetic 
(H) fields can be decomposed into TE and TM polar- 
izations, described by scalar fields tp satisfying Neumann 
(TE, i' = H^) and Dirichlet (TM^ = E;,) boundary 
conditions at the metallic surfaces |13 |. 

To analyze these geometries, we employ two comple- 
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FIG. 1: Schematic of the cyhnder-plate geometry consisting 
of two infinite metallic cylinders of radius R, separated by a 
distance a from each other, and at a distance h from one (top) 
or two (bottom) metallic sidewall(s). 



mcntary and exact computational methods, based on 
path integrals (PI) and the mean stress tensor (ST). The 
methods are exact as they involve no uncontrolled ap- 
proximations and can yield arbitrary accuracy given suf- 
ficient computational resources. They are complemen- 
tary in that they have different strengths and weaknesses. 
The PI method is most informative at large separations 
where it leads to analytical asymptotic expressions. The 
ST method, while relatively inefficient for large separa- 
tions or for the specific geometries where PI is exponen- 
tially accurate, is formulated in a generic fashion that al- 
lows it to handle arbitrary complex shapes and materials 
without modification. As both of methods are described 
in detail elsewhere [isl . we only summarize them 
briefly here. The present geometry provides an arena 
where both methods can be applied and compared. 

In the PI approach, the Casimir force is calculated via 
the constrained partition function. The Dirichlet (Neu- 
mann) constraints on the TM (TE) fields are imposed by 
auxiliary fields on the metallic surfaces [l5l | which can be 
interpreted as fiuctuating charges (currents). The inter- 
action between these charges is related to the free-space 
Green's functions — the addition of metallic sidewall(s) 
merely requires using image charges (currents) to en- 
force the appropriate boundary conditions. The calcula- 
tions are further simplified by using Euclidean path inte- 
grals and the corresponding imaginary-frequency uj = iw 
Green's function. In the case of infinite cylinders, these 
surface fields can be represented in terms of a spectral 
basis: their Fourier series, leading to Bessel functions in 
the Green's function [1, [13, [l3| • An important property 
of such a spectral basis is that its errors go to zero expo- 
nentially with the number of degrees of freedom. 

We also use a method based on integration of the 
mean ST, evaluated in terms of the imaginary-frequency 
Green's function via the fiuctuation-dissipation theo- 
rem [Tsj. The Green's function can be evaluated by a 
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FIG. 2: Casimir force per unit length between two cylin- 
ders (black) vs. the ratio of sidewall separation to cylin- 
der radius h/R, at fixed a/R = 2, normalized by the to- 
tal PFA force per unit le ngth b etween two isolated cylinders 
[FpFA = |(fic7rVl920)i/R/V [lil]. The solid lines refer to 
the case with one sidewall, while dashed lines depict the re- 
sults for two sidewalls, as indicated by the inset. Also shown 
are the individual TE (red) and TM (blue) forces. 



variety of techniques [13| , but here we use a simple finite- 
difference frequency-domain method (FDFD) that 
has the advantage of being very general and simple to 
implement at the expense of efficiency — it is much less ef- 
ficient for this specific geometry than the PI method. The 
results from both methods are shown in Fig. [21 with the 
PI method indicated by solid/dashed lines for two/one 
sidewalls, and the ST method indicated by data points. 
Both results agree to within the numerical accuracy, as 
expected, although we have fewer data points (and larger 
error bars) from the ST method because it is less efficient 
for this geometry. We focus on the interpretation of the 
results rather than on the computational techniques. 

To begin with, we compute the force between the two 
cylinders (with a/R = 2) as a function of the sidewall 
separation h/a, for fixed a. The results, for both one 
sidewall (solid lines) and two sidewalls (dashed lines) are 
shown in Fig. [5] for the total force (black), and also the 
forces for the individual polarizations TE (red) and TM 
(blue). The forces are all normalized to the PFA result 
between two cylinders , which is independent of h and 
does not affect the shape of the curves in Fig. [21 Two 
interesting observations can be made from this figure. 
First, the total force (for both one and two sidewalls) is 
a nonmonotonic function of h/R: at first decreasing and 
then increasing towards the asymptotic limit between two 
isolated cylinders. The extremum for the one-sidewall 
case occurs at h/R~ 0.27, and for the two-sidewall case 
is at h/R Ki 0.46 (similar to the h/R Ki 0.5 for squares 
between two sidewalls [^). Second, the total force for 
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FIG. 3: Casimir force per unit length between two cylinders of 
fixed radius R vs. the ratio of sidewall separation to cylinder 
radius h/R (for one plate), normalized by the total PFA force 
per unit length between two isolated cylinders. The force is 
plotted for different cylinder separations of a/R = 0.2, 0.6, 
1.0, 2.0, 3.0, and 4.0. 



the two-sidewall case in the h = hmit is larger than 
for /i — > oo. As might be expected, the /i-dcpendcnce 
for one sidewall is weaker than for two sidewalls, and the 
effects of the two sidewalls are not additive: not only is 
the difference from h ^ oo force not doubled for two 
sidewalls compared to one, but the two curves actually 
intersect at one point. 

Since nonmonotonic sidewall effects appear to occur for 
a variety of shapes (both for square Q and circular cross 
sections) , it is natural to seek a simple generic argument 
to explain this phenomenon. As we see in Fig.[2l and also 
in our earlier work 0], the nonmonotonicity arises from 
a competition between the TE and TM force contribu- 
tions: the TE force is quickly decreasing with h while the 
TM force is slowly increasing. Therefore, an underlying 
question is why the TE force increases as the sidewall 
comes closer, while the TM force decreases. 

An intuitive perspective of the effects of the metallic 
sidewall(s) on the TE/TM forces is obtained from the 
"method of images," whereby the boundary conditions at 
the plate(s) are enforced by appropriate image sources. 
For the Dirichlet boundary conditions (TM polarization) 
the image charges have opposite signs, and the potential 
due to a charge (more precisely, the Green's function at 
any imaginary frequency, which determines the Casimir 
force) is obtained by subtracting the contribution from 
the opposing image. Any configuration of fluctuating TM 
charges on one cylinder is thus screened by images, more 
so as h is decreased, reducing the force on the fluctuating 
charges of the second cylinder [2^ . Since the reduction 
in force is present for every configuration, it is there also 
for the average over all configurations, accounting for the 
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FIG. 4: Casimir force per unit length between a plate and 
two cylinders of fixed radius R vs. the ratio of cylinder 
separation to cylinder radius a/R, normalized by the total 
PFA force per unit length between a cylinder and a plate 
FppA = |(fic7rV960)\/E/2a7 The force is plotted for 

different plate separations of h/R = 0.28, 0.6, 1.0, and 2.0. 
Note that the normafization is different from the cylinder- 
cylinder PFA in the previous figures. 



variations of the TM curves in Fig. [21 By contrast, the 
Neumann boundary conditions (TE polarization) require 
image sources (current loops) of the same sign. The total 
force between fluctuating currents on the cylinders is now 
larger and increases as the plate separation h is reduced. 
(An analogous additive effect occurs for the classical force 
between current loops near a conducting plane.) 

Note, however, that while for each fluctuating source 
configuration, the effect of images is additive, this is 
no longer the case for the average over all configura- 
tions. More precisely, the effect of an image source on 
the Green's function is not additive because of feedback 
effects: the image currents change the surface current 
distribution, which changes the image, and so forth. For 
example, the net effect of the plate on the Casimir TE 
force is not to double the force as ^ 0. The increase is 
in fact larger than two due to the correlated fluctuations. 

While the method of images explains the competition 
between TE and TM modes that underlies the nonmono- 
tonic effects, further considerations are required to ensure 
that their sum is nonmonotonic. For example, if the TE 
and TM variations with h were equal and opposite, they 
would cancel with no net dependence on h. That this 
is not the case can be checked by examining the limit 
a ^ h ^ R. In this limit the forces arc dominated by 
the lowest spectral (Fourier) mode, s-wave for TM and 
p-wave for TE pTj]. The former is stronger and leads to 
an asymptotic form (for one plate): 
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confirming the reduced net force as the cylinders ap- 
proach the plate. While the logarithmic dependence on 
R could have been anticipated the h'^ scaling is a 
remarkable consequence of the multi-body effect. Each 
cylinder and its mirror image can be considered as a 
dipole of size ~ h. The interaction of the two dipoles 
should scale as the interaction between two cylinders of 
size h with Neumann boundary conditions. For a ^ h 
the force for the latter prob lem scales as ^ h'^/a^, ex- 
plaining the above result [lOl , up to the logarithm. To an- 
alytically establish the nonmonotonic character, we also 
need to show that the TE force is dominant in the op- 
posite limit of h R. So far, we only have numerical 
arguments in favor of this . 

In Fig. El we show the total force vs. h/R for a variety 
of different values of the cylinder separation a in the pres- 
ence of a single sidewall. The value a/R = 2 corresponds 
to our previous results in Fig. [21 Note that if a is too 
large or too small, the degree of non-monotonicity (de- 
fined as the difference between the minimum force and 
the h = force) decreases. (For small a, the force is 
accurately described by PFA, while for large a the TM 
mode dominates as indicated in Eq. H]).) The separation 
a/R = 2 from Fig. |2l seems to achieve the largest value 
of non-monotonicity. 

When the force between the cylinders is not monotonic 
in h, it also follows that the force between the cylinders 
and the sidewalls is not monotonic in a. A nonmonotonic 
force Fx between the cylinders means that there is a value 
of h where dF^/dh = 0. Since the force is the derivative 
of the energy, F^ ~ ~d£/da, at this point d'^S/dadh ~ 0. 
These two derivatives, of course, can be interchanged to 
obtain d{d£/dh)/da = 0. But this means that dFy/da — 
at the same point, where Fy = —d£/dh is the force 
between the cylinders and the sidewall. This cylinders- 
sidewall force is plotted in Fig. |4l as a function of a/R 



for various values of h/R. (It is not surprising that the 
effect of a small cylinder on the force between two bodies 
is smaller than the effect of an infinite plate. This is also 
reflected in the fact that the cylinder-cylinder force is 
generally less than the cylinder-plate force for the same 
cylinder diameter and surface separation [Tsj.) 

The advantage of the cylinder-plate force compared 
to the cylinder-cylinder force is that it seems opera- 
tionally closer to the sphere-plate geometries that have 
been realized experimentally. In order to measure the 
cylinder-cylinder force, one would need to suspend two 
long cylinders in vacuum nearly parallel to one another. 
To measure the cylinder-plate force, the cylinders need 
not be separated by vacuum — we expect that a similar 
phenomenon will arise if the cylinders are separated by a 
thin dielectric spacer layer. Unfortunately, the nonmono- 
tonic effect in Fig. [His rather small (roughly 0.2%), but it 
may be possible to increase this by further optimization 
of the geometry. In future calculations, we also hope to 
determine whether the same phenomenon occurs for two 
spheres next to a metal plate. 

In previous research, unusual Casimir force phenom- 
ena were sought by considering parallel plates with ex- 
otic materials: for example, repulsive forces were pre- 
dicted using magnetic conductors [l^, combinations of 
different dielectrics [l^l, fluids between the plates [2l[, 
and even negative- index media with gain [22l |. A dif- 
ferent approach is to use ordinary materials with more 
complicated geometries: as illustrated in this and previ- 
ous [9| work, surprising nonmonotonic (attractive) effects 
can arise by considering as few as three objects. 
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